Quantum quenches and work distributions in ultra-low-density systems 
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We present results on quantum quenches in systems with a fixed number of particles in a large region. We 
show that the typical differences between local and global quenches present in systems with regular thermo- 
dynamic limit are lacking in this low-density limit. In particular, we show that in this limit local quenches 
may not lead to equilibration to the new ground state, and that global quenches can have power-law work dis- 
tributions ("edge singularities") typically associated with local quenches for finite-density systems. We also 
show that this regime allows for large edge singularity exponents beyond that allowed by the constraints of the 
usual thermodynamic limit. This large-exponent singularity has observable consequences in the time evolu- 
tion, leading to a distinct intermediate power-law regime in time. We demonstrate these results first using local 
quantum quenches in a low-density Kondo-like system, and additionally through global and local quenches in 
Bose-Hubbard, Aubry-Andre, and hard-core boson systems in the low-density regime. 



Introduction. Motivated by remarkable experimental 
progress in realizing and exploring non-equilibrium physics 
in cold-atom systems [1|, there has been increasing interest 
in the dynamics of thermally isolated systems . Despite 
the rapidly growing body of research in this class of non- 
equilibrium dynamics, many aspects are still poorly under- 
stood. For example, what type of equilibration can be ex- 
pected for various types of local and global quenches? An- 
other question involves the overlap distribution, closely re- 
lated to the work distribution 13-5] for a quantum quench. 
What is the typical form of the distribution of overlaps of the 
initial state with the final eigenstates? What are the effects 
of various overlap distributions and work distributions on dy- 
namical (time-evolving) quantities? 

In the experimental settings suitable for exploring non- 
equilibrium physics, such as cold atoms and semiconductor 
nanostructures like quantum wells, a common situation is to 
have a fixed number of particles in a large spatial region. This 
contrasts sharply with the solid-state notion of the thermody- 
namic limit, where large regions are filled with a constant den- 
sity. The study of non-equilibrium issues (e.g., quenches and 
work distributions) in such situations, where the usual ther- 
modynamic limit is not applicable, is clearly of topical impor- 
tance but has been near-absent in the non-equilibrium theory 
literature. 

In this work, we focus on this ultra-low-density limit — 
fixed number of particles, arbitrary large sizes. We present a 
study of quenches in a system which is the counterpart of the 
Kondo model in this low-density regime. We present several 
dynamical aspects which, through calculations in a few other 
low-density systems, we show to be generic features of quan- 
tum quenches in this limit. 

One peculiarity of this limit is a blurring of differences 
between the consequences of local versus global quenches. 
Another striking result involves the overlap distribution, 



is the initial state 



(tf(0)| 4>)l } ) , where |*(0)) = 
(ground state of initial Hamiltonian), and m indexes the eigen- 
states of the final Hamiltonian. We show that this quan- 
tity is dominated by a power-law decay, ^mT a , genetically 
for quenches involving low-density systems. The associated 
"edge singularity" in the work distribution has large power- 
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FIG. 1. Kondo-like model with single mobile fermion. (a) Energy 
gap as a function of J, L — 100 sites, showing distinct behaviors in 
three regimes. Boundaries between A,B (B,C) regimes are denoted 
according to criteria f = L (£ = 1). (f = localization length.) 
Left inset shows system geometry. Right inset shows spectrum in C 
regime, (b) Density profiles rij in three regimes; dashed lines show- 
ing exponential localization. 



law exponents which would not be compatible with the usual 
thermodynamic limit. This in turn has remarkable conse- 
quences on the real-time evolution: in the evolution of ob- 
servables away from their initial value, there appears an in- 
termediate power-law regime between the initial perturbative 
time period and the large-time steady-state behavior. 

Since the fixed-number large-size limit is applicable to 
many experimental non-equilibrium setups, these results are 
expected to be relevant to experimental situations realized or 
realizable in the near future. 

Kondo-like model. The main system we use for demon- 
strating these general results involves a few (N c ) mobile 
fermions ("conduction electrons") in a tight-binding closed 
chain (Fig. fTh. One site of the lattice is Kondo-coupled to a 
single spin-^ "impurity". The Hamiltonian is 



H 



c\ s Cj +M + h.c. 



J Simp ' "0 (1) 



where S = Y^ s 



&ss' cq,s' is the s P m on site i = (s, s' 



are spin indices), and i £ [0, L — 1] is the site index. We study 
quenches of J, i.e., local quenches, starting from the ground 
state at J = Jj and studying the dynamics after changing J 
instantaneously to its new value Jf. The ground state is a 
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FIG. 2. Local quench in Kondo-like model with N c = 1 mobile 
fermion. (a) Time evolution of occupancy no at impurity-coupled 
site; C->C quench, Ji = 10 2 , J f — 10; L — 100 sites. The long- 
time average {ho)DE (solid black line), around which no(t) oscil- 
lates, is significantly different from the equilibrium J = Jf value, 
(b) (ho)DE as function of Ji, for final values in A,B,C regimes: 
Jf = 0.02 (green), 0.9 (dark green), 10 (blue). Dots mark the pa- 
rameters for panel c. In (a,b), gray dashed and red dash-dotted lines 
show equilibrium no values for J = Ji and J = Jf. (c) Size depen- 
dence of (no)r,E values for Ji — 10 2 . Upper blue: Jf — 10; lower 
green: Jf = 0.9. Red dots are corresponding J = Jf equilibrium 
values. 
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FIG. 3. Overlap distribution 

a, and work distribution p(oj) 
(a,b) Quenches from Ji — 10" 



) , power-law exponents 
(a-d) Kondo-like model, N c = 1. 



^1 



to Jf = 10 3 



J/ = 10 3 (C^C), Ji = 10 3 to Jf = 10~ 3 (C 



(A->C), Ji = 10 to 
-+A), and Ji = 10~ 3 



to Jf = 10~ (A— >A). (c,d) Exponents of power-law fits m~ a 
for quenches ending at and starting from the A regime. (e,f) Bose- 
Hubbard chain with Nt = 3 bosons in L — 20 sites; interaction 
U quenched from 0.3 to 0.5. (f) Approximations to the work distri- 
bution are obtained using Gaussians of width a to replace the delta 
function of Eq. ([2j. (g) Kondo-like model, N c — 3; quench from 
Ji = 10" 3 to J f = 10 _1 . 



spin singlet, and quenches of J preserve the spin, so that all 
dynamics is confined to the spin singlet sector. 

In Fig. [JJ we summarize the equilibrium physics of the 
N c = 1 system. In an infinite chain, in the ground state, the 
fermion is localized around the impurity-coupled site (i = 0) 
with localization length £. (£ decreases with increasing J.) At 
large J (regime C), the itinerant fermion is almost completely 
localized at site (£ < 1). At smaller J, the itinerant fermion 
is spread over multiple sites £ > 1 (regime B). In an infi- 
nite system, this region would extend to arbitrarily small J. 
However, for any finite size L, there is a boundary-sensitive 
small- J regime (regime A) where the fermion cloud extends 
over the whole system (£ > L). For 1 < N c <$lL, we naturally 
get additional features, but the same general behavior persists 
in the three regimes. 

Observables. We will present time dependences of the oc- 
cupancy nrj(t) of site i = for the Kondo-like system, and 
of the Loschmidt echo C{t) = \ (*(0) |*(t) ) | 2 . The observ- 
able no(t) is of obvious importance for the model ([TJ, while 
C(t) is well-defined for any model and is closely related to 
the work distribution [3]-[5]. Despite the nonlocal nature of 
the Loschmidt echo, there exist proposals for experimentally 
measuring this quantity, and related quantities have been mea- 
sured J6)- 

Lack of equilibration to new ground state in local quenches. 
The final value at which an observable O saturates is given by 



df) 



O 



the so-called 



(O)de = £ m <*(0) 
"diagonal ensemble" (DE) value Q. In Fig. |2|a) we show 
the time dependence of n after a quench within the C region. 
We note that reaches the DE value (uo)de relatively 

rapidly, and then shows 'revivals' at roughly periodic intervals 



of t ~ L/2, The DE value where n saturates is markedly 
different from the ground state value of no for J = Jf. This 
seemingly contradicts the intuition that a local quench in a 
large system should lead to relaxation to the final ground state 
value, because the energy pumped into the system by a local 
quench is a effect. The reason this does not happen in 

the C— >C quenches is that the itinerant electron only occupies 
a small number of sites near the impurity position. Thus, most 
of the lattice sites do not play any role in the dynamics, and 
cannot serve as a bath to absorb the disturbance at site 0. This 
effect is not restricted to N c = 1, but is true for finite number 
N c > 1 of fermions for L — > oo ||8l. 

In Fig. |2|b) we show (?t.o)de as a function of Ji for fixed 
Jf. The (?t.o)de values deviate significantly from the J = Jt 
equilibrium values for most Jj, Jf combinations. Fig. |2jc) 
shows, through L-dependences of (tio}de> tnat the lack of 
equilibration in quenches to C or B regions is not a finite-L 
effect. 

This effect represents a loss of the distinction between local 
and global quenches, which is a generic feature of the L — > oo 
limit with finite particle number. 

Overlap Distributions. Fig. |3ja-d) summarize overlap 
distribution behaviors in quantum quenches between different 
regimes of the system ([TJ for N c = 1. These behaviors can be 
derived from detailed consideration of the eigenfunctions 

In C^C quenches, the ground state overlap (<f>Q | </>|£L 

much larger than the others, while the small m ^ overlaps 
have the form oc sin(cim) 0. The most remarkable feature is 

the power-law behavior, (<j) ^ | ^JP) ~ mT a ', in quenches 

starting from or ending in the A region. The exponent a is 2 
for A^A quenches and 1 for A— s-C quenches. 
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These power law behaviors are a generic phenomenon; we 
have found such power-law overlap distributions in several 
other systems in the low-density limit, both for local and 
global quenches. (The behavior is particularly clean for the 
N c = 1 system because of its simplicity.) Fig. [3jg) shows 
the overlap distribution for the same model with N c = 3 
fermions. There are now additional structures, but the dom- 
inant overlaps follow a clear power law. Fig. [3|e) shows the 
overlap distribution for a Bose-Hubbard chain at low density 
0. Again, there are interesting additional structures, but the 
dominant overlaps follow a clear power law (~ mr a ~). 

Work distribution. The overlap distribution is related to 



(2) 



where 



E. 



(/) 



( f) 

Eq are the final eigenenergies mea- 
sured from the final ground state energy. This is the so-called 
work distribution GH5], except for a shift between cj and the 
usual work variable. (The energy prior to the quench plays 
no role in the temporal dynamics and so is not relevant for 
this work.) The work distribution is related to the Loschmidt 
echo: C(t) = \J dujp(uj)e lu,t \ 2 . Since£(0) = 1 by definition, 
p{oj) must be normalizable. 

At large sizes (but constant particle number), p(co) can be 
treated as a continuous function starting from uj = A, the 
finite-size gap, which vanishes at large L. We have found that, 
for quantum quenches in low density systems, the work dis- 
tribution generically has behavior p (lj ) « Pau~ b for lu > A, 
with large exponents b > 1. These power-law divergences are 
analogs of what would be called "X-ray edge singularities" in 
systems with a regular thermodynamic limit. In finite-number 
systems, p(co) remains normalizable despite the singularity as 
A — > because the magnitude of p(oj) also vanishes (po — » 0) 
in the large-size limit, due to the vanishing density. 

This contrasts sharply to systems with the usual thermo- 
dynamic limit where density remains constant as L — >• oo, 
and p(lo) itself is a well-defined non-vanishing quantity in the 
limit. This constrains the singularity p(oS) ~ to have 
smaller exponent, b < 1 (e.g., J4)). The low-density systems 
of interest here have no such constraint; a central result of the 
present work is that super-linear singularities (b > 1) are sig- 
natures of low-density systems. 

For the model (fTJ with N c = 1, p(u) — cj~ 5/2 (A-sA) 
and p{w) ~ w _3/ MA^C). (b > 1 in both cases.) Fig. ^f) 
shows the work distribution for a global interaction quench in 
the Bose-Hubbard chain, with the delta function regularized 
as gaussian. There is a power law with super-linear {b > 1) 
singularity. This is another example of the loss of distinction 
between global and local quenches in the low-density limit, as 
"edge singularities" are normally associated only with local 
quenches for finite-density systems [4|. 

We have also found super-linear singularity exponents 
in other low-density systems [8(, e.g., quenches of the 
strength/position of a weak trapping potential for a Bose- 
Hubbard system, quenches of on-site potentials and hopping 
strengths for hard-core bosons in a ladder geometry, and 
quenches of quasi-disorder potential strengths in an Aubry- 
Andre [9| system. 
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FIG. 4. (a-f) Time evolution of no(t) and C(t) in iV c = 1 model, 
(a-d) A->A quenches: Jt = IO" 3 , J f = 1(T 2 . (e.f) A->C: J l = 
10~ 3 , Jf = 10. Extended intermediate region between ^t 2 region 
and long-time oscillatory region is seen in no(t) for A—)- A quenches 
and in C(t) for both A— >A and A— »C quenches, but not in no(t) 
for A— >C quenches, (g) Bose-Hubbard, global interaction quench, 
3 bosons in chains of length L — 10 and L — 30. Quench from 
Ui = 0.02L to U f = 0.01L. (f) N c = 3 fermions; J; = 10~ 3 , 
J f = 10" 2 . 



Role of the density of states. For the model ([T| with N c — 
1, the behavior ^mT a implies energy-dependence ^oj~ a < 2 
for the overlap distribution. Together with a factor of w -1 / 2 
from the ID single-particle density of states, this leads to 
p(u) - oj-*- 1 / 2 , i.e., b = 5/2(3/2) for A-^A(C) quenches. 

This argument can be generalized: if the overlap distri- 
bution follows vnT a and the density of states in the rele- 
vant lower-energy part of the spectrum behaves as p(u) ~ 



the work distribution p(uj) 



will have exponent 



b = 2ja — 7 + 2a (8J. For single-particle systems, we have 
7 = —1/2 in ID, as in the above example. For a generic sys- 
tem, however, the many-body density of states does not nec- 
essarily behave as a power law. We have found cases (Bose- 
Hubbard chain with trap) where an approximate power-law 
region with exponent 7 in p(uj) leads to an approximate power 
law in p(uj) with exponent b = 2^a — 7 + 2a (8). Also, if 
a = 1/2, any power-law form of p(uj) implies a linear edge 
singularity p(iS) ~ In this case, a super-linear edge 

singularity can only happen with some non-power-law form 
of p{oj). This occurs in the Bose-Hubbard chain case of Fig. 

ie,f) m. 

The intermediate-time ~ t@ region. The appearance of 
larger powers in the edge singularity has novel consequences 
for real-time dynamics. We have identified an intermediate- 
time power-law region in the dynamics of the Loschmidt echo 
(and other observables), that appears as a direct consequence 
of the large-power edge singularity. 

At initial times after a quench, observables and C(t) evolve 
away from their initial value quadratically with time, ~ t 2 , 
as can be explained from generic perturbative arguments. We 
have found that, when p(u) has a large-exponent singularity, 
there is a region of time (after the initial perturbative times 
and before the large-time steady-state oscillations), where 
£(0) — £(t) = 1 — C(t) follows a new power-law behavior. If 



p(uj) 



' with b G (1, 3) in the energy range [A, A] and the 
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contributions outside this energy window can be neglected, 
then in the time window between t ~ A -1 and t ~ A -1 one 
sees the behavior 1 - £(t) ~ i 6 " 1 0. 

The same phenomenon is also found in some observables: 
\0(t) — O(0)| can have an extended region after the initial 
perturbative ~ t 2 region with a new exponent ~ t@o . When 
has the form of a rank-1 projector, = \x)(x\> we can 
write the time evolution as 0(t) = J efcjp^a^e*" , where 

P< 5( W ) = £ % - I ^X^ 5 I X) (3) 

m 

differs from the work distribution pi in that one factor of 

( f) — 

the overlap is replaced by (<j>m \ X)- If Po( w ) h as a power- 
law singularity structure u>~ b ° with exponent 6^ € (1,3), 
the time evolution of 0(t) away from 0(0) will show the 
intermediate-time region ~ t^o (J}q = — 1). When the 

operator does not have the form O = |x)(x|, it is not sim- 
ple to formulate an analogous expression. A generic operator 
for a many-body system will not have this form, but the site 
occupancies for systems with single itinerant particles (e.g., 
rij for our N c = 1 system) have the forms of rank-1 projec- 
tors, as does the Loschmidt echo for any system. Currently, 
little is known about Pq(ui) behaviors for projector- type ob- 
servables in different quantum quenches, or about the con- 
ditions necessary for having an intermediate-time regime in 
generic observables. 

The intermediate-time regime for low-density systems is il- 
lustrated in Fig. |4] For the model ([T]i with N c = 1 parti- 
cle, this regime is present in the Loschmidt echo, for both 
A^A and A— >C quenches. In the occupancy n (t), the 
intermediate-exponent regime can be seen for A— > A quenches 
(with form ^t 1 / 2 ), but not for the A— >C quenches, for which 
case the eigenstate dependence of (q^iP \ x) does not favor 
a large enough exponent in pn (uj) [8|. Fig. |4^g) displays 
the intermediate-time region for a Bose-Hubbard chain with 
interaction quenches, and Fig. |4|h) shows the same for the 
Kondo-like model with N c = 3 itinerant fermions. 

With hard-core bosons on a ladder-shaped lattice, consid- 
ering time evolution after various local and global quenches 
of hopping strengths and on-site potentials, we find extended 
intermediate-time regions ~ t p in 1 - C(t), with exponents 
matching j3 = b — 1 where b is the singularity exponent 
in p(u), calculated with gaussian regularization (8). With 
quenches of a trapping potential, we find quench parameter 
combinations where p{uS) shows super-linear edge singular- 
ities (b > 1) but no intermediate-time regime shows up in 
the £(t) dynamics because the singularity exponents are too 



large, b > 3 J8J. We have also found an example (Aubry- 
Andre system) where there are well-defined p(us) ~ oj~ b re- 
gions but the contributions from outside the power-law region 
are so large that the dynamical intermediate-time signature is 
washed out (8). 

Extent of the intermediate-time region. If the power-law 
window for p(oj) is u £ [A, A], the ^i* 3 region with j3 £ 
(0,2) extends from t ~ A -1 to t ~ A -1 . The scale A is 
generally of the order of the bandwidth, and so is set by the 
hopping strength. Since the finite-size gap A vanishes with 
increasing system size L, the intermediate-time region gets 
more and more extended in time for larger L. This is shown 
in Fig.Qg) through a comparison of two different L values. 

Discussion. For systems that are not well-described by 
the traditional thermodynamic limit but instead have a fixed 
number of particles in a large size, as is common in setups 
relevant for large classes of non-equilibrium experiments, we 
have presented universal features of quantum quenches. These 
include edge singularities with large exponents not possible 
in 'regular-limit' systems, a loss of the usual distinctions be- 
tween local and global quenches, and a novel intermediate- 
time region in the dynamics. 

Universal behaviors in quantum quenches are generally 
sought and discussed in asymptotic times. A new universality 
at intermediate times, visible in widely different systems, is of 
obvious distinction and interest. 

Our results open up various questions and research avenues. 
One issue is to bridge the gap between the regime considered 
here and the regular thermodynamic limit. A systematic study 
of overlap/work distributions and associated quench dynam- 
ics with varying density and system sizes is currently lacking. 
A related issue is that of experimental accessibility. Since a 
measurement requires finite density, it is important to demar- 
cate which densities in a real finite-size system show features 
of our fixed-number large-size limit, from those correspond- 
ing to the regular thermodynamic limit. It may also be inter- 
esting to supplement our results on the Loschmidt echo with 
time-evolution studies of traditionally measurable observables 
such as densities and correlation functions. Finally, Eq. ([3]) 
highlights the general lack of knowledge about expectation 
values of observables in (and overlaps of states with) gen- 
eral eigenstates of many-body Hamiltonians. Current activ- 
ity on the Eigenstate Thermalization Hypothesis is addressing 
eigenstate expectation values of some observables |7][T0|, but 
clearly further investigations are warranted. 

Acknowledgments. MH thanks M. Vojta for discussions 
on equilibrium properties of the Hamiltonian ([TJ. 
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Supplementary Materials 



S.I. OVERVIEW 

In these Supplementary Materials, 

• We provide a derivation of the existence of an extended 
intermediate-time power-law region in the Loschmidt 
echo C(t), when the work distribution p(uj) has a 
power-law behavior uj~ b with exponent b G (1,3). 
(Section[SJll) 

• Considering the case where the density of states has 
power-law behavior, we derive a relation between vari- 
ous exponents. (Section S.III ) 



We give some details for quantum quenches in four 
low-density systems. (Interaction quench in Bose Hub- 
bard chain, trap quench in Bose Hubbard chain, sev- 
eral quenches in hard-core bosons on a ladder geom- 
etry, and quasiperiodic potential quench in an Aubry- 
Andre lattice.) We show various examples of super- 
linear p(ui) ~ co~ b behaviors, and examples of the 
intermediate-time power-law region in 1 — £(t). (Sec- 
tion |SlVl ) 

We show time evolution for the Kondo-like model with 
N c = 3 fermions, to show that there is no equilibration 



S.V (Only N c = 1 is shown in the main text.) 



i We provide additional details and derivations for local 
quenches in the Kondo-like model for N c = 1, (Section 

CT ) 



S.II. DERIVATION OF THE INTERMEDIATE REGIME 

FOR£(i) 

In this Section, we prove that an edge singularity p{uS) ~ 
u)~ b with b £ (1,3) leads to the intermediate time behavior 
^t b_1 for £(0) — C(t), provided that contributions of parts of 
p(ui) outside this power-law region can be neglected. 

We consider p(io) to have power-law behavior p(u) ~ i>J~ h 
in the energy window ui G [A, A]. Here A is an energy 
scale of the order of the finite-size gap, which should van- 
ish in the i^oo limit. The energy A up to which the power 
law holds is system-dependent, generally of the order of the 
bandwidth. There are two time scales A -1 and A -1 corre- 
sponding to these energy scales. We show below that, when 
the exponent b is in the appropriate range, b £ (1,3), the 
novel intermediate-time regime 1 — C(t) ~ t h ~ x appears in 
the time window between these two timescales, when contri- 
butions from u ^ [A, A] are neglected. 

The Loschmidt echo C(t) is related to the Fourier transform 
of the work distribution: 

G(t) = [ p{uj)e itu duj (S.l) 



as£(i) = G*(t)G(t). 
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Neglecting the contributions from energies outside the 
power-law region, we obtain 



G(t) oc / u- b e ltuJ duj 

J A 

= A 1 - h £ b (-itA)-A 1 - b S b (-itA), (S.2) 

where £b(z) = e~ ZUJ uj~ b duj is the generalized exponen- 
tial integral (6, z £ C). When 1 — C(t) is plotted using 
the above expression ( |S.2) , we see a clear intermediate-time 
^t b ~ 1 region between a ^i 2 and a ~t°, as shown in Fig. SI 
for b = 2.5. 




FIG. SI. Time evolution of th e Loschmidt echo C(t) = G(t)G*(t) 
obtained from expression i S.2 1 with b = 2.5, A = 10~ 4 , A = 4. 



Below we derive the behaviors of C(t) at short and interme- 
diate times using expansions of the exponential integral £b(z). 



S.II.-l. Short time regime. 



When tA <C 1, both terms in ( |S.2| i can be expanded in 
Taylor series: 



a £b{—itA) = {-ity-'rii -b) + 



A 



1-6 



6-1 



A 2-b A 3-b 

; 6-2 <_ 2(6-3)' 



+ «i j^-ttti ^^ 2 + 0(t 3 ), (S.3) 



and 



A^^i-itA) = (-tty-'Til - 6) + 



b-ii 



A 



1-6 



6-1 



A 



2-6 



A 3- 



+ 2(6^3) <2 + ° (t3) - (S4) 



Using (|S.3|>,([ST4[> in ( |S.2| i, and assuming 6 e (1,3) and 
A -1 < A -1 , we get 



£(t) = |G(i)G*(t)| = 1 + |^—^A h - 1 A 3 - b t 2 



C(i 3 ). (S.5) 



Here G(t) was normalized using the criterion £(0) = 1, 



S.1I.-2. Intermediate-time regime. 

In the intermediate time interval A -1 <C t <C A -1 the 
expression for A 1-6 £f,(— it A) is still given by the Taylor ex- 
pansion ( |S.3| l, but we cannot use ( |S.4| i because tA is not small. 
Instead we use the large-time asymptotics for £ b (— it A): 



A L -"£ b (-itA) = e 



HA 



iA- 



Oit^A- 1 



(S.6) 



Case 1: 1 < 6 < 3. When 6 is in this range, using the ex- 
pression \S3\ for small tA and (|S.6|l for large tA, we obtain, 



.1-6 



G(t) oc - T(l- 6) cos [|(1- 6) 



,6-1 



,2-6 



6-2 



t - «r(l - 6) sin 



£<!-»> 



6-1 



(S.7) 



As before, G(t) has to be normalized such that £(0) = 1. 
Then the Loschmidt echo at leading order is 



C(t) = l+2(l-6)r(l-6) cos ^-b) 



A 6 " 1 * 6 " 1 . (S.8) 



Case 2: 6 > 3. In this case the leading term is (tA) 2 
rather than (tA) b_1 : 



G(t) oc 



A 



1-6 



.3-6 



6-1 2(6-3) 



- t 6 " 1 cos 



.2-6 



|(1 - 6)] T(l- 6) -i^t. (S.9) 



This gives, for the Loschmidt echo at leading order, 

C(t) = 1 b — 1 A 2 * 2 

() (6-3)(6-2) 2 ' 



(S.10) 



Therefore there is no observable intermediate-time region 
with exponent different from 2. 



S.III. RELATIONSHIP BETWEEN POWER-LAW 
EXPONENTS IN OVERLAP DISTRIBUTION, WORK 
DISTRIBUTION, AND DENSITY OF STATES 

In cases where the density of states has a power-law behav- 
ior, p(uS) ~ cl> 7 , the eigenstate index scales with the energy 
(relative to ground state) as 



dJp(wf) ~ w 7+1 . 



(S.H) 



Therefore for an overlap distribution with m a behavior we 
will have the dependence 



1 4P) 



-2a( 7 +l) 



(S.12) 



The definition of the work distribution p(oS) contains a factor 
of the density of states in addition to this overlap squared; thus 
the exponent 6 of the edge singularity p(cj) ~ cj~ b is 



6 = 27a + 2a 



1- 



(S.13) 
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In the case of the N c = 1 Kondo-like model in the main text, 
we had 7 = — 1/2 (single-particle density of states in one 
dimension), and a = 2 (1) for an A to A (C) quench. This 
leads to super-linear singularities in the work distribution: b = 
a + i = 5/2 (3/2) for the two quench regimes. 

A consequence of ( |S.13| l is that, if the exponent for the over- 
lap happens to be a = 1/2, the work distribution p(uS) will 
have linear singularity (b = 1) for any power-law form of the 
density of states. 

Note that, in a generic multi-particle situation, we have no 
a priori reason to expect a clean power-law behavior of p(cu) 
at small u. In Section [5. 1 V.A| we will see a case where, even 
though a = 1/2, it can combine with a p(uS) with some more 
complicated behavior to give ap(w) with super-linear singu- 
larity (b > 1). 



S.IV. VARIOUS MODELS: 
SUPER-LINEAR EDGE SINGULARITIES 
AND 

INTERMEDIATE-TIME REGIME 

In this section, we present numerical results on local and 
global quenches in several different systems in the low-density 
limit. This sampling of low-density systems shows that 
power-law behaviors in the overlap distribution, and super- 
linear power-law behaviors in the work distribution [p(u>) ~ 
oj~ h with b > 1], are generic in this important limit. The 
signature in real-time dynamics (an extended intermediate- 
time power-law regime), requires the additional condition b < 
3 and that contributions outside a single power-law regime 
in p(ui) can be neglected. Therefore, the intermediate-time 
regime appears in some but not all cases. 



In S.IV. A we describe results for a global quench of the 
interaction strength, for a Bose-Hubbard chain in the dilute 



limit. In S.IV.B we consider Bose-Hubbard systems in a weak 
harmonic trap, and consider quenches of the trap strength and 
the trap position. In |S.IV.C| we present results for a dilute 
system of hard-core bosons in a ladder geometry, considering 



both local and global quenches. Finally, in S.IV.D we consider 



a single particle in a quasi-periodic (Aubry-Andre) potential, 
and present results for quenches of the strength of the quasi- 



periodic ('disorder') potential. For the systems of S.IV.C and 



S.IV.D we also present some examination of the effect of in- 
creasing system size. 



S.IV.A. Bose Hubbard, interaction quench 



We consider the celebrated Bose-Hubbard Hamiltonian 

L-l L 

^bh = - J2 6 5 ft i+i + 5>ifo - x ) (S - 14) 

in one dimension. Here bj, b\ are bosonic operators for site j, 



and n,j — bjbj are site occupancies. Energies (times) are mea- 
sured in units of the hopping (inverse hopping) strength. We 



use open boundary conditions and consider global quenches 
of the interaction parameter U, form Ui to Uf. 

Figure S2 shows numerical results for Nf, = 3 bosons in 
L = 30 sites, which we expect to represent the "fixed Nf, in 
large L" limit. We show results for initial and final values 
of U in both the small (U < 3) and large (U > 3) regimes. 
The phenomena of large-exponent power-law singularity in 
the work distribution and intermediate region in the time evo- 
lution of C(t) appear for both large and small quenches, and 
for final parameter in either the large-C/ or the small-[/ regime. 



The overlap magnitude 



dfh 



plotted in the top 



panels, naturally shows more scatter and more structures com- 
pared to the N c — 1 Kondo-like model treated in the main 
text. However, in each case there is a dominant set of data 
points which behave as ~ m -1 / 2 . In the large U cases, where 
there are several 'bands', there are also large contributions 
from higher bands. The extra features (the large-energy con- 
tributions, and the contribution of the many low-energy states 
for which the overlap is nonzero but falls below the dominant 
to -1 / 2 states) will of course leave signatures in the tempo- 
ral dynamics. However, our signature phenomenon (extended 
intermediate regime between perturbative and final regimes) 
appears in all these cases. 

In the second and third rows, we show continuous approx- 
imations for the density of states p(w) for the final U = Uf, 
and the work distribution p(ui). This is done by replacing the 
delta functions by gaussians of energy width a, 



8{u 







1 



27TCT 



exp 



2<T 



in the definitions 



p(u) =^5(uj - e m ) 

m 



(S.15) 

(S.16) 
(S.17) 



Here e m — Em — Eq 1 are the final eigenenergies measured 
from the final ground state energy. 

we show continuous curves with a = L _1 for 



?(f) 



In Figure 



S2 



p(u) and with u = 3L _1 for p(oj). The choice of a is a com- 
promise for visualization; with smaller a one sees more oscil- 
latory behavior associated with the discreteness of the spec- 
trum, while choosing larger a washes out features near the 
beginning of the spectrum. 

The density of states p(ui) has a maximum in the central 
region of each band. For present purposes, the relevant part 
of the spectrum is where p(uj) is increasing toward its first 
maximum; this is naturally highlighted in the log-log plots of 
the second row. 

In the smoothed work distributions (third row), the region 
uj < er shows a broad plateau which is an artifact of our 
smoothing procedure, and looks artificially broad on a loga- 
rithmic scale. We have mostly omitted this part from the plot- 
range. After this part (u > a), there is an extended spectral 
region where p(ui) follows an approximate power law. The 
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FIG. S2. Nt = 3 bosons in L — 30 sites. The interaction U is quenched; the (Ui,U{) pair is indicated on top of each column. Top row: 



overlap distributions. The dashed line in second panel is 



-1/2. 



the dominant overlaps follow same exponent in every case. For large 



Uf (middle two panels) significant contributions are visible from higher bands. Second row: many-body density of states p(lj) at U = Uf . 
Insets show same with frequencies in linear scale. At large Uf, there are several bands. Third row: Work distributions p(oj). Dashed lines 
are ~b~ 3 ^ 2 and ~6 -2 . Approximate power-law behavior of p(w) clearly extends into ui values near the top of the (first) band, where p{uj) 
definitely cannot be approximated with a power-law. For large U / (middle two panels), the contributions from higher bands are not visible in 
the energy range shown. Bottom row: Evolution of Loschmidt echo away from its initial value. A region with intermediate exponent between 
the initial ~t 2 and the final oscillatory is visible in all cases. 



exponent b of this "edge singularity" p(u>) ~ u b is super- 
linear, with b m 2 for small Uf and b slightly smaller for large 
U f . 

The approximate power-law behavior of p(u>) emerges in 
a more complicated manner than in the N c = 1 system of 
the main text. In this case, if the relevant part of the den- 
sity of states followed a power law p(uS) ~ w 7 , it would not 
be possible to have a super-linear edge singularity because 
a = 1/2 would imply b = 1 for any 7 (see Section S.IIIi. 
Remarkably, the non-powerlaw form of p{uS) conspires with 
the m~ 1 / 2 overlap distribution to cause an approximate power 
law in p(oo), with exponent b > 1. At present it is not clear 
whether the power-law behavior of p(lo), or the exponent val- 
ues, becomes exact in some limit. 



The p(uj) 



behavior results in the dynamical feature 



that we have highlighted in this work: in the bottom row, we 
see clear intermediate regimes between the perturbative and 
the steady-state regimes. The behavior in this regime is ap- 
proximately ~ t° , with j3 w 1 for small Uf and j3 slightly 
smaller for large Uf. This is consistent with our prediction 
of = b — 1 (Section S.IIi. Note that, in the large Uf cases, 
there are oscillations in the ~ t* 9 region, which result from the 



contribution of the higher bands. 

It would be interesting to find out if an extended interme- 
diate power-law region is visible in the time evolution of ob- 
servables that are more commonly studied (or more likely to 
be experimentally measurable) than the Loschmidt echo. 

As the density (filling) is increased, the features we have 
presented gradually disappear. Fig. [S3] shows the case of unit 
filling (8 bosons in 8 sites). The density of states and work dis- 
tributions are significantly different; we do not analyze them 
here. 

Curiously, in weak quenches between small values of U 
(leftmost panels), there is a sequence of dominant over- 
laps that seems to follow a mT 1 ! 2 behavior, and a small 
intermediate-time region does appear in 1 — C(t). In the cor- 
responding work distribution (not shown), there is no obvi- 
ous p(uj) ~ u)~ b behavior. It is somewhat surprising that the 
intermediate-time regime appears at unit filling. To explain 
this, one could consider two different types of L — > 00 limits. 
Presumably, the small intermediate-time region will disappear 
in the large-i limit taken with filling held constant (usual ther- 
modynamic limit), but would get more extended in the large- 
L limit taken with fixed Nf, = 8. A systematic study of the 
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FIG. S3. Nh — 8 bosons in L — 8 sites. The interaction U is quenched; the (Ui,Uf) pair is indicated on top of each column. Top row: overlap 
distributions. The dashed line in first panel is ~m _1//2 . There seems to be a power-law behavior in this weak quench, but not in any of the 
other quenches. Bottom row: Loschmidt echo. A small intermediate-exponent region is seen in the weak-quench between small U values, but 
not in any of the other quenches. 



shape of the overlap distribution, work distribution and den- 
sity of states as a function of filling and size is clearly called 
for, but is beyond the scope of this work. Some L-dependence 



is explored in Sections S.IV.C and S.IV.D 



No intermediate-time regime is seen in the dynamics for 



the other quenches in Fig. S3 nor is there any power law like 
behavior in the overlap distributions. This is expected because 
the super-linear edge-singularities and the intermediate-time 
regime in C(t) are novel features of low-density systems. 



S.IV.B. Bose Hubbard in harmonic trap; trap quench 

We now consider bosons on an open-boundary chain sub- 
ject to to a harmonic confining trap in addition to the Bose- 
Hubbard Hamiltonian ( |S.14| >: 



H 



BH+trap 



Hi 



BH 



(S.18) 



Harmonic traps are fundamental to considerations of cold- 
atom experiments, the most prominent experimental setting 
for non-equilibrium dynamics in isolated systems. 

We consider quenches of both the trap strength fc tr and the 
trap center jo- The Hamiltonian ( |S.18| l has various regimes of 
possible interest, such as large and small interaction U, strong 
and weak trapping potential, etc. For present purposes, we 
will confine ourselves to small U and weak trapping, k tr ~ 
0(1O -3 ). The weak and strong trapping regimes are loosely 
analogous to the A and C regimes of the Kondo-like model 
detailed in the main text. 

In Figure S4 we show data for a trap strength quench 



(^tr 



0.001 to fc L 



(/) 



strength quench (j'q 
the overlap distribution 



= 0.0012) and for a trap position 
7.3 to j ( Q f) = 8.3). In both cases 
' 1 1 shows a dominant se- 
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FIG. S4. Bose-Hubbard in trap: trap strength quench (left) and 
trap position quench (right). Left: Nt = 4 bosons in L = 19 sites; 
interaction U = 0.2; trap center at jo = 7.3; trap strength quenched 
from k$ = 0.001 to fc^ = 0.0012. Right: N b = 3, L = 28; 
U — 0.2; fc tr = 0.002. Trap center quenched from jg = 7.3 
to j — 8.3. Dashed lines in top row (overlap panels): ~m ' 
and ~m -2 . Dashed lines in second row (density of states panels): 
~ cj. Dashed lines in third row (p(co) panels): ~oj -5 and ~oj~ 7 . No 
intermediate-exponent region is visible in the lowest row (Loschmidt 
echo evolution). 



a is between 3/2 or 2. As in the case without a trap (Section 
S.IV.Ai, we show approximations to the density of states p(oj) 



ries which follows an approximate power law; the exponent 



obtained by replacing delta functions with gaussians. 

The density of states is similar to the case without a trap 



10 



(S.IV.Ai, but the small-w behavior of p(uS) seems closer to a 



power-law form in the presence of a trap. This exponent (7) 
is close to 1. The argument of S.III then predicts the work 
distribution p(cu) ~ us~ b with 6 = 5 (for a = 3/2) or b = 7 
(for a — 2). Indeed, p(oj) does have an approximate power- 
law decrease with exponent around the range 5 to 7, in the 
energy window where p(ui) increases roughly linearly. It is 
difficult to be more certain about the exact values of the ex- 
ponents, because of the scatter in the overlap distribution and 
because of the uncertainties associated with replacing delta 
functions by gaussians in order to plot p(oj) and p(ui). Also, 
it is difficult to be certain from the present data whether there 
is a genuine power-law region in p(uj) or whether it is merely 
approximate. Despite the difficulties with rigorous determina- 
tion of exponents, the relation ( |S.13| l is at least approximately 
valid. 

In this low-density system and these types of quenches, we 
thus have edge singularities in p(u>) with exponents far larger 
than that allowed in systems having he usual thermodynamic 
limit. However, as explained in Section S.II an exponent in 
p(u>) ~ Lu~ h with b > 3 does not lead to a distinguishable 
intermediate-time region in the evolution of the Loschmidt 
echo. Indeed the bottom panels of Figure S4 show the ^t 2 



region directly followed by the oscillatory region. The singu- 
larity exponent is too large to see an intermediate-time sig- 
nature in C(t). Of course, an extended intermediate region 
might be visible in some observable other than the Loschmidt 
echo. An exploration of the evolution of physical observables 
in such quenches remains an open task for future work. 



was obtained by Gaussian regularization of the delta function 
(as in Sections [STVAl and |SiVB) i with width a = 0.1. For 
the system sizes considered here, p(cu) is noisy near the band 
edges as there are not enough states in this regions to obtain a 
smooth curve. The same phenomenon appears in the work dis- 
tribution p (w) shown in Fig. S5 d). This quantity is sharply 
peaked around the bottom of the spectrum. The power-law 



lines in Fig. S5 d) are only guides to the eye. 

Even though it is difficult to assign unambiguous power- 
law exponents forp(w), the existence of an intermediate-time 



regime in the temporal dynamics is clear, as seen in Fig . S5 e) 



S5 



d), are 



for 1 — C(t). The exponents in p(uj) ~ w , Fig 
consistent with the exponents for the interme diate -region ( 
t 13 ), Fig.[S5|;e), following the results of Section|S.Il| /3 = b- 1. 



'e), following the results of Section ! 
Fig. S5 e) also shows the perturbative region for small t and 
the oscillatory behavior obtained at large times. The time at 
which the large-time os cillato ry behavior starts (correspond- 
ing to ~ A -1 of Section S.II I clearly grows with system size. 
The extent of the intermediate-time regime gets progressively 
wider as the system size increases at fixed particle number. 



Finally Fig. S5 g), shows how C(t) itself (not subtracted 
from £(0) = 1) looks like in the intermediate-time region. 



S.IV.D. Single particle in Aubry-Andre potential 

We consider a single particle on an i-site chain subject to 
a quasiperiodic potential, i.e., a cosine potential with incom- 
mensurate period: 



S.IV.C. Hard-core bosons in a ladder 

We next consider a system of Nb hard-core bosons on a 
ladder of N v vertices with vertex-dependent on-site potentials 
Vi and bond-dependent hopping Jij : 



L-2 



HCB 



E J a h t h i 



E^ 



(S.19) 



where b\ and bi are hard-core bosonic operators for site i with 



h , b 



1 and b\ 2 — b\ = 0. One leg of the ladder is taken 

to have one site more than the other, in order to avoid spurious 
symmetries. 

In the initial Hamiltonian, the hopping terms are set to 
J.- L j = — 1 and the on-site potentials are set to Vi = . At 
t = some of the hopping terms or on-site potentials are 
changed abruptly according to the color coding of Fig. \S5[a) 
where Vi — is colored gray and = — 1 yellow. We con- 
sider a local quench of a rung hopping strength (Jy) (first col- 



umn of Fig. S5 I, a global quench of the rung hoppings (second 
column), and a local quench of an on-site potential Vi (third 
column^ 

Fig. 



S5 



as a function 



b) shows the overlaps 

of the eigenenergies e m . The overlaps have a maximum mag- 
nitude at 772 = for all the quenches considered. The density 



of states p(u) of the final Hamiltonian, given in Fig. S5 c), 



H = -J 



E 



C ] + 1 C 3 



E^ 



(S.20) 



with Vj = V cos (27rgij) having irrational wave vector qi, 



here taken to be q\ 



VE-i 



This is known as the Aubry- 



Andre model. There is a localization transition at V/J = 2; 
single-particle eigenstates are exponentially localized in space 
for V/J > 2. We will consider quenches within the delocal- 
ized regime, from V = 0.1J to V = 0.6 J. 

As V is changed from toward 2 J, the tight-binding band 
splits up into sub-bands and continues splitting further until 
it becomes fractal at V — 2 J. Fig. S6 a) shows the density 
of states for V — 0.6J, where the (sub)band edges are seen 
through \u> — w c dgc|~ 1//2 cusps characteristic of single-particle 
bands in ID. (The energy gaps between (sub)bands are some- 



times too small to be seen.) Fig. S6 b) shows the work dis- 
tribution p(uj). The work distribution shows edge singulari- 
ties at every band edge and is thus quite intricate. Fig S6 c) 
shows that the work distribution near the bottom of the first 
band (w— >0) has a clear power law behavior; p (us) oc w~ 3 / 2 . 
As expected, in order to fulfill the normalization condition 
J dujp (u>) = 1, the magnitude of p(u>) decreases with L. The 

d)]. There 



S6 



pre-factor po decreases as L~ 4 overall [Fig. 
are some fluctuations because of the quasi-random nature of 
the system. The gap b etwe en the ground and the first excited 
state, A, shown in Fig. S6 e), shows a characteristics L~ 2 de- 
cay with L. 
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FIG. S5. Three different quenches for a system of Nb — 3 hard-core bosons in a ladder geometry with different ladder sizes, (a) Geometry 
of the system and quench type. Yellow bonds correspond to Ji.j — — 1 and gray vertices to Vi = 0. Blue bonds and black vertices show 
parameters changed in the quench, (b) Eigenstate overlap as function of energy relative to final ground state energy, (c) Density of states 
normalized to the total number of states, (d) Work distribution. Black dashed lines are ~ uj~ b with exponents b chosen by eye. (e-f) 
Loschmidt echo. In (e), the perturbative (~ t 2 ) and intermediate ~ t 6-1 regions are highlighted with dashed black power-law lines. 



Because p(io) has large contributions in several regions out- 
side the small-cj power-law region, no intermediate-time re- 
gion is visible in the time evolution of the Loschmidt echo. 



S.V. LACK OF EQUILIBRATION IN KONDO-LIKE 
MODEL WITH N c = 3 FERMIONS 



time evolution of the local occupancy at the impurity-coupled 
site, nrj(i), after a quench from J, = 10 to Jf = 100. The 
long-time average (no)DE is the black solid horizontal line 
around which no(t) oscillates. This is markedly different from 
the final equilibrium value, shown as the red dot-dashed hori- 
zontal line. 



In the main text we presented the absence of equilibration 
to the final equilibrium value for N c = 1 fermion, in contrast 
to the expectation for local quenches in finite-density systems. 
This is not a single-particle curiosity, but a generic feature of 



low-density systems. In Fig. S7 we show an example of the 



S.VI. DETAILS FOR iV c = 1 KONDO-LIKE MODEL 

In this Section we provide details for results presented in 
the main text for the Kondo-like model with a single mobile 
fermion (N c = 1). In IS. VI. A| we provide derivations of the 
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FIG. S6. Single particle in Aubry- Andre potential, (a) Density of states p{ui). (b) Work distribution p [uf] for different system sizes L. 
(c) Rescaled work distributions at small u. Values of po are chosen such that pg 1 p(uj) for different L lies on the same curve. A power-law 



dependence, p(oS) 



with b — 3/2, is very clear, (d) Scaling of the pre-factorpo as a function of L. (e) Dependence of the gap A on L. 




the singlet sector: 



li> = 4(cSjn>|t)-c} t |fi)|;>) (S.21) 



V2 



where j is the site index and |0) is the fermionic vacuum. 



S.VI.A. Derivation of the overlap behaviors 



FIG. S7. Blue solid curve: Time evolution of the local density no(t) 
in the Kondo-like model, with N c = 3 fermions, L — 20, J; = 10, 
Jf = 100. Solid black line: the long-time average value of no(t), 
i.e., the diagonal average value {nojcs- The gray dashed and red 
dot-dashed lines indicate the ground-state values of no correspond- 
ing to J = Ji and to J = J/. 



forms of the overlap distributions presented in Figure 3 of the 
main text. In S. VLB we show how parameter values for A^A 



We summarize below analytic calculations for the overlaps 



quenches can be chosen in a way such that power law region of 
the overlap distribution gets systematically extended with sys- 
tem size. In S.VI.C we explain why for A— >C quenches there 



is no intermediate-time region visible in nn(t), even though 
such a region is visible in C(t), as seen in Figure 4 of the 
main text. 

Spatial index as basis label. Since we are restricted to the 
singlet sector, for N c — 1, we can use the position of the mo- 
bile fermion as the label for a complete set of states spanning 



for quench cases starting and ending in A and 



(/)\ 

m / 

C regimes. 

We start with A— ^A quenches, J, t <C 1. We derive the 
overlap behavior (power law, ^m~ 2 ) by treating Jj t pertur- 
batively. Next, we derive the overlap behaviors for quenches 
between A and C regimes when J;(j)<Cl and Jfuy^l, by us- 
ing zeroth-order expressions for the eigenfunctions in the two 
limits. For C— >C quenches, we derive the sine behavior of the 
overlap by using the zeroth-order ( J3>1) expression for the 
final eigenstates. 

In the main text, the index m was used to label the eigen- 
states within the symmetry sector where the z-component of 
the total spin is zero. Within this symmetry sector, there are 
also triplet states and states with odd spatial reflection parity, 
which have zero overlap with the initial state. Here, we re- 
strict further to eigenstates with nonzero overlap, which are 
about one-fourth of the states spanned by the m index. We 
use pL to index this restricted set of eigenstates. 
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FIG. S8. Overlap distributions for the N c = 1 Kondo-like model 
with L — 100, compared with analytic predictions. Only the even- 
parity excited singlet states are indexed by fi. (a) A— >A quench, 
Ji = 10" 3 , Jf = 10"\ Red solid line is Eq. j§j25). (b) A ^C 
quench, Ji = 10" 3 , Jf = 10 3 . Red solid line is Eq. ( (S^g). (e ) 
C^A quench, Ji = 10 3 , Jf = 10" 3 . Red solid line is Eq. dS.29t . 
(d) C^C quench, Ji = 10 2 , Jf = 10 4 . Red solid line is Eq. l |S3Tt . 



S.VI.A-1. A^-A quench 



This expression is compared to exact numerical values in Fig. 
la). 



S. VI.A-2. A->C quench 

The final ground state is not included in the power-law 
regime for A to C quenches, reflecting the fact that the ground 
state in the C regime is well-separated and quite different 
from the other eigenstates. In the calculations below, we are 
therefore only interested in \i > 0. For very large J (ex- 
treme C regime), the particle is localized at site i = in the 
ground state, and excluded from this site in the other eigen- 
states. Thus the excited eigenstates ("band states") are the 
single-particle wave functions in a chain of length L — l with 
hard- wall boundary conditions: 



,c,(o) , 
Wo * 



L-l 

E^ 7 

i=i 



sin 



L 



(2fj, + 1) j \j). (S.26) 



The initial ( J <C 1) state is 
ire tl 

V2 



,A(0) 



(0)) =i" 1 / 2 E,k)at 



lowest order. Therefore the overlaps are 

L-l 



E sir 



2?r 



(S.27) 



For Jj ,/<SCl, we write both initial and final Hamiltonians as 
Hi j = H o + JijV and use first-order perturbation theory. 

The reflection-symmetric eigenfunctions of the unper- 
turbed Hamiltonian Hq are the plane waves: 



L-l 



E^ 1 



i=0 



and the energies are 
eigenfunctions are 



(0) 



# (1) > 



L/2 

# (0) > + E 



-2 W -J |j) (S.22) 
2 cos [j/i] . At first order, the 



4 A,(0) 



s (o) _ E (o) 



^ (0) ; 



(S.23) 



In the limit L — > oo one can replace summation with integra- 
tion, leading to 



C (0) l«o 0) ) 



2^2 



(S.28) 



tt(2/z + 1) 

This expression is compared to exact numerical values in 



Fig. S8 ;b). 



5. VZA-3. C—>A quench 



For ^ 3> 1 and J/ <C 1, we use the zeroth-order expres- 
sions in the two limits: q£ ) — 4>q'^ ) = J2i=o ^*,o \ x ) 



and 



Therefore, the overlaps of the ground state 



qb^ ) of Hi with overlap distribution is the constant function 



i f) : 



given by Eq. ( |S.22| i. The resulting 



the excited states 



) of Hf are, at first order, 



£' (0) 



^ (0) > 



(S.29) 



-V, 



(o) 



J/ — J 

^(0) ^(0) V,o 



(S.24) 



This expression is compared to exact numerical values in Fig. 

ISSJc). 



V 



are found to be 



The matrix elements l/^o ■ 

constant, V^.q = 3\/2 (4L) . Thus the overlaps as function 
of n are 



3^2 |J/-Jj| 
81 1- cos[2ttL-V] 



(S.25) 



S.VI.A-4. C-^C quench 

We prove the sine behavior for quenches from anywhere in 
B or C regimes to the extreme C case (Jf — oo). Describing 
the sine behavior more generally for anywhere in C to any- 
where in C involves the same physical ideas, but is too clumsy 
in notation to write out. 
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For the final excited states we use the zeroth order expres- 
sion ( |S.26| l. For the initial ground state, we can formally 



use the exact form 



b { l) ) = A/"" 1 / 2 cxp [-[x]/£(Ji)], where 

[x] = min(x, L — x) is the distance from the impurity coupled 
site. Here £(J) is a decreasing function of J and M(J) is a 
normalization factor. This gives for the overlaps 



4" 



b [f) ) 



2V2 



L/2 



sm 



(2/1+1) 3! 



(S.30) 

-|x| 



IS 



By definition of the C regime, £(Jj) < 1 and exp 

a fast decaying function. Hence the main contribution to the 
sum is given by the first term (x = 1): 



In Figure |S9| we plot A— sA overlap distributions for sev- 
eral system sizes. As L is varied, we adjust J,- t such that the 
ratio of £(J) to L stays fixed. Here £(J) is the localization 
length that the mobile fermion would have in case of an infi- 
nite system. In the A regime, £(J) > L. Adjusting Jij in 
this manner ensures that we do not get into the B regime as 
L is increased. In addition, it turns out that increasing L and 
decreasing Jj j in this coordinated manner yields a series of 
curves that systematically extends the power-law behavior, as 
seen in Figure [S9{a). 



Figure 

6$° I 4® 



S9 



shows the logarithmic derivative of the 

versus m data, which extracts the power-law 

exponent. We note that the region where this quantity is nu- 
merically a = 2 (dashed horizontal line), also gets systemati- 
cally extended through this procedure. 



2V2 



(2/i + 1) 



(S.31) 

This expression is compared to exact numerical values in Fig. 
IHb). 



S. VLB. Scaling demonstration of the power law in A to A 
quench 



S.VI.C. Absence of intermediate-time regime in rij (t) for A to 
C quench 

For A— >C quenches, an extended intermediate-exponent 
regime in the time evolution is seen in the Loschmidt echo 
£(£) but not in the local density no(i). To clarify this situ- 
ation, we note that hj is a rank-1 projection operator: hj — 
Using the language introduced in the main text for 
such observables, instead of the work distribution we should 
consider 

p A . H = % - E m )(<t>$> I 1 3) (S.32) 




10 1 10 2 
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FIG. S9. The overlap distribution \ 



) and the numerically 



determined power-law exponent (absolute value of the logarithmic 
derivative), plotted for different sizes. Sizes L = 100, 200, 400, 800 
are used with the Jij adjusted so as to keep £,{Ji,f)/L fixed. For 
L = 100 we used J- = 10 -3 , J f = 10~ 2 . 



The A^A quench overlap data in Figure 3 of the main text 



4? 



dfh 



for 



clearly shows a power-law behavior 

a particular system size L. However, since the definition of the 
A region is itself i-dependent, the extension of this behavior 
with increasing L is not immediately obvious. 



One factor of the overlap distribution is now replaced by 
(4>m^ I j)- I n tne C region, for the excited ('band') eigenstates 
rn > 0, this quantity is approximately constant except when 
it vanishes for symmetry reasons. In the excited eigenstates 
m > 0, the quantity is very small for the impurity-coupled 
site j = 0, and 0(1/ L) for the other (j > 0) sites. In either 
case, this leads to the power-law behavior pf lj (u>) cx uj~ hfl j 
with 6ft. = 1. This exponent is not large enough to cause a 
distinct intermediate-time region in \rij(t) — rij(0)\. 

In contrast, the exponent relevant for the Loschmidt echo 
is b = 3/2, appearing in the work distribution. This leads to 
an intermediate-exponent region 1 — C(t) ~ t 1 / 2 between the 
initial ^t 2 and the large-time oscillatory regimes. (Figure 
4 in main text.) 

Note that, in addition to having an unsuitable exponent, the 
quantity pn (w) corresponding to site j = has large contri- 
bution from m = 0, which is outside the power-law region. 
The corresponding quantities for other sites, say pa L/2 (w) for 
j = L/2, would have larger spectral weight in the power-law 
region. However, since 6 rlj = 1 for all sites j, there is still no 
intermediate-time region for j ^ 0. 



